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Abstract 

We analyze the determination of volume effects for correlation functions that depend on an ex- 
ternal momentum. As a specific example, we consider finite volume nucleon current correlators, 
and focus on the nucleon magnetic moment. Because the multipole decomposition relies on 5*0(3) 
rotational invariance, the structure of such finite volume corrections is unrelated to infinite volume 
multipole form factors. One can deduce volume corrections to the magnetic moment only when 
a zero-mode photon coupling vanishes, as occurs at next-to-leading order in heavy baryon chiral 
perturbation theory. To deduce such finite volume corrections, however, one must assume continu- 
ous momentum transfer. In practice, volume corrections with momentum transfer dependence are 
required to address the extraction of the magnetic moment, or other observables that arise in mo- 
mentum dependent correlation functions. Additionally we shed some light on a puzzle concerning 
differences in lattice form factor data at equal values of momentum transfer squared. 
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I. INTRODUCTION 



QCD in the non-perturbative regime has proven notoriously difficult to understand quan- 
titatively. This non-Abelian gauge theory of gluons coupled to quarks is simple enough to 
write down, and even becomes perturbative at high energies. At low energies, however, 
the theory is strongly coupled which results in the confinement of quarks and gluons into 
hadrons. Accounting for the measured properties of hadrons is a challenge; making reliable 
predictions is even more difficult. After more than a quarter century of dedicated work, lat- 
tice gauge theory has made considerable progress in addressing strong interaction physics. 
As a first principles numerical technique, lattice gauge theory can be used to determine low- 
energy hadronic properties rigorously from QCD. For a comprehensive overview of lattice 
methods, see (l| 

Despite considerable success, lattice QCD calculations still suffer from a number of sys- 
tematic errors. Due to computational restrictions, simulations cannot be carried out at 
the physical pion mass. Instead heavier pions must be used. Additionally current lattice 
volumes are not significantly larger than the typical hadronic length scale, which for most 
observables is set by the pion Compton wavelength. To eliminate such systematic errors, 
it is desirable to have an independent tool that predicts the pion mass, and lattice volume 
dependence of observables. Fortunately, effective field theory techniques exist suited for 
this purpose. Chiral perturbation theory, for example, provides a model independent and 
systematic tool to address the pion mass and lattice volume dependence of many low-energy 
hadronic properties. Pioneering work on the volume dependence of the chiral condensate 
and pion observables appeared a while ago in series of papers @, MMM- Additionally finite 
volume amplitudes have been analyzed in order to study unstable particles and multiparticle 
physics j6|, 3, HI- A variety of observables have since been treated in finite volume. Recent 
work has included treating heavy mesons and baryons P, [HI QQE3, extending the 
validity of finite volume corrections for meson masses and decay constants [3, EEI, la] , ex- 
ploring different regimes in finite volume theories [l7, 18, 19], and investigating multiparticle 



systems in finite volume [20|, |2l|, |22j, |23 



turn 



Further work atfinite volume has included processes depending on an external momen- 

volume corrections to electromagnetic observables, for ex- 
29| . A point that was specifically addressed in these latter 



25, 26, 27 



Deducing 

ample, can be a subtle task [28j 
works was the lacking connection of finite volume amplitudes to low-energy multipole expan- 
sions. Multipole expansions are inherent in the description of SO(3) invariant physics. At 
finite volume such a description of low-energy matrix elements ceases to be valid. Utilizing 
effective field theory, finite volume corrections to amplitudes can be computed, but not cor- 
rections to multipole moments and other observables, such as multipole polarizabilities. Said 
another way, an effective field theory only allows one to match the calculation of correlation 
functions. Thus chiral perturbation theory allows one to calculate the model- independent, 
long-distance behavior of finite volume correlation functions. Connection to infinite volume 
observables, which might be accessible through various different correlation functions, is not 
mandated. 

In this work, we consider finite volume correlation functions depending on an external 
momentum. Specifically, we detail the finite volume modifications to nucleon current matrix 
elements, and analyze whether it is possible to deduce volume corrections to the nucleon 
magnetic moment. Based on symmetry considerations, we show that such a connection to 
infinite volume physics is not generally possible; but, within the framework of heavy baryon 
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chiral perturbation theory, the vanishing of a new finite volume coefficient at next-to-leading 
order allows for the connection to be made. We are careful to expose, however, that this 
analysis relies on the assumption of continuous momentum transfer. On the lattice, the 
available momentum modes are quantized, but there is a regime in which the assumption of 
continuous momentum is approximately valid. This regime, however, is beyond the reach of 
current computing resources. Thus in practice, the volume effect must be determined from 
the momentum transfer dependent nucleon current matrix element. 1 Although we have spe- 
cialized to current matrix elements, our analysis generalizes to any momentum-dependent 
correlation function at finite volume. Finally using our expressions for current matrix ele- 
ments, we show that form factors calculated at equal values of momentum transfer squared 
(but at differing momentum transfer) can differ due to volume effects. This potentially 
resolves a puzzle seen in lattice form factor data. 

Our presentation has the following organization. In Sect. [Til we write down the single 
particle effective action for the nucleon coupled to an electromagnetic field. This is done 
for a spatial torus. Here we show based on symmetry arguments that the electromagnetic 
current is additively renormalized, and that there is an additional magnetic-like zero-mode 
interaction. Next in Sect. IHH we calculate polarized nucleon electromagnetic current matrix 
elements in finite volume. Heavy baryon chiral perturbation theory is utilized, and the 
relevant finite volume functions appearing in the calculation are listed in Appendix [A] The 
results of this calculation can be used to determine the volume- dependent coupling constants 
in the single nucleon effective action. In Sect. IIVI we investigate the conditions necessary 
for the matching to be performed, and find they require lattice sizes larger than currently 
available. In Appendix [Bj we carry out the momentum expansion one non-trivial order 
further to discuss the magnetic radius at finite volume. The form factor difference puzzle is 
taken up in Sect. HVCl Generic features of our findings are summarized in Sect. [V] 



II. SINGLE PARTICLE EFFECTIVE ACTION 



We begin by considering the low-energy dynamics of a nucleon in an external electro- 
magnetic field. We use a two component isospinor N for the nucleon. In infinite volume the 
nucleon can be described by the Lagrangian 

C = N (iv ■ D - - ^g- \S„, SJ F>A N, (1) 
V 2M N 2M N J K ' 

where v 11 is the heavy nucleon four- velocity, is the covariant spin operator, and D± = D^— 
v^{y ■ D), see j33|. The coefficient of the kinetic term is exactly fixed by reparametrization 
invariance j34[, but the form of the operator is unique only up to field redefinitions. The 
gauge covariant derivative appearing above is = + ieQ^A^, where Qn is the nucleon 
charge matrix. The coupling constant fi^ is a diagonal matrix containing the proton and 
neutron magnetic moments. In writing the above Lagrangian, we have included all zeroth 
and first order terms in an expansion in photon frequency. At higher order in the frequency 



1 We mention that background field methods alternatively utilize two-point functions to deduce magnetic 



30, 



31 



3211. As we focus 



moments from the shift in particle energies linear in the magnetic field, see e.g. 
on three-point functions, our analysis has no direct relation to finite volume corrections in background 
field calculations. 
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expansion, there are terms with more powers of the field strength tensor, F^ v ', and derivatives 
thereof. Because we work well below the pion production threshold, we have integrated out 
pion-nucleon interactions to arrive at Eq. ([I]). The low-energy constants in the matrix /xjv 
depend on the pion mass and other couplings of the nucleon theory with pion interactions. 
That theory is chiral perturbation theory, and we treat such dependence as implicit. 

In finite volume, we can write down a theory analogous to Eq. ([T]). To be concrete, we 
assume that this theory is defined on a torus of length L in each of the three spatial directions. 
This reflects the underlying lattice QCD action with quarks subject to periodic boundary 
conditions, but ignoring any possible effect from the finite time direction. The general 
form of the single nucleon effective action can be constructed by writing down all possible 
operators consistent with the symmetries of electromagnetism on a torus. The discrete 
symmetries C, P, and T remain. Boost invariance and 5*0(3) rotational invariance possessed 
by the Lagrangian in Eq. (jTJ), however, are reduced to the cubic symmetry group, which is 
isomorphic to 5*4. On compactified spaces, gauge transformations are more constrained 
for the gauge field zero mode. A well known example of this occurs in finite temperature 
gauge theory, see, e.g. (35j. On the spatial torus we consider, the zero modes of the three- 
vector potential have periodicity constraints under gauge transformations. This restriction 
on gauge transformations leads to the ability to construct more gauge invariant operators 
than in infinite volume. 



In [28], the spatial analog of the Polyakov line was used to write down gauge invariant 
operators involving the gauge field zero modes. Specifically employed was the Wilson line 
Wj defined by 

Wi = expfj jdxiAA , (2) 

where the line integral starts at some point x^ and runs to x^ + Lxi. There is no implied sum 
over % on the right-hand side of Eq. (j2J). Because the gauge potential is periodic, 2 cycling 
the compact dimension must produce a gauge invariant object. The factor of e/3 reflects 
that the quark charges are quantized in such units. Under a gauge transformation, we have 
— > + d^a. In order for the Wilson line in Eq. ([2]) to be gauge invariant, we must have 

67T 

a{Xi = L) = a{Xi = 0) H n i; (3) 

e 

where is an integer. As this condition must hold in each spatial direction, we find 

a(x) = —jj 1 ' x + a(x), (4) 

where a(x) is a periodic function of the The linear term in the gauge function a(x) gives 
rise to a quantized shift of the gauge field zero mode. 

To write an effective theory of the nucleon and photons in finite volume, it is easiest to 
use linear combinations of Wilson lines that are Hermitian and have definite C, P, and T 



2 In lattice simulations that probe electromagnetic observables, either a background electromagnetic held is 
gauged into the action, or current matrix elements are calculated. In the latter case, the current operator 
J^{x) = -0(2;) 7^(2;) is periodic. Hence in an effective theory for such calculations, one must take to 
be periodic otherwise the effective action is not single valued. Our discussion, however, does not apply to 
the background field method which requires a separate treatment altogether. 
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transformations. These are the operators and w\ ^ given by 



= ~ (w t + Wj) , (5) 
Wi _) = i fWi - Wj) • (6) 



2i 

It is easy to show that these two operators are not independent, and we choose to work 
with the because they have the same C, P, and T symmetry transformations as 

the Ai. Constructing the most general gauge invariant Lagrangian on a torus is arbitrarily 
complicated. As our interest is with the nucleon magnetic moment, we restrict our attention 
to operators having only one insertion of \ Multiple insertions of Wj" ^ lead to photon- 
nucleon couplings with more than one photon. Furthermore we restrict our attention to 
operators with at most one derivative. 3 We will investigate the conditions that justify such 
a photon frequency expansion in Sect. HVl 

We now build the single particle effective action for nucleons and photons on a torus. To 
write this theory, we abandon the covariant notation employed in Eq. (P) above, because the 
torus has only S± invariance. Including all single Wf~' ) operators with at most one derivative, 
we have the Lagrangian 

D 2 CAL) - , v a N (L) C 2 (L) ( r A 

2M N 2M N 2M N 2M N \ J 



£ = N 



N. 

(7) 

where iD = i (^D — D^j defines an Hermitian operator. Notice that VjWf "* = 0. In writing 

Eq. (J7|), we have explicitly indicated the dependence on the spatial length L. The coupling 
constants above run with the infrared cutoff, l/L, and must be determined by matching 
calculations in the microscopic theory. For example, the magnetic moment operator is 
accompanied by the finite volume coefficient Hn{L) which is given by 

H N {L) = n N + Sfx N (L), (8) 

where the term S^n{L) is the finite volume effect that can be determined using chiral per- 
turbation theory. Running the infrared cutoff to zero produces the infinite volume magnetic 
moment, lim Hn(L) = fM^. Compared to Eq. (jTJ, there are two new operators allowed 

L— >oo 

by symmetry. These operators contain single photon couplings as well as a tower of cubic 
invariant multi-photon couplings. There are, however, further multi-photon operators that 
we have not written in Eq. (JTJ). Such operators involve multiple insertions of \ The 
new coupling constants C\(L) and C^L) appearing in Eq. ([7j) both must run to zero when 
the infinite volume limit is taken. As we work below any multi-particle thresholds, these 
L-dependent couplings run to zero exponentially fast in asymptotically large volumes 
Notice that expanding the Lagrangian in Eq. (J7]) to linear order in the gauge field, we have 
an accidental 5*0(3) invariance. 

The effective theory in Eq. (0) can be used to calculate single photon-nucleon processes. 
For example, for an unpolarized nucleon at rest, we have the current matrix element 

(N(0)\J»\N(0)) = Q N eg»°, (9) 



3 Terms with further derivatives are considered in Appendix [B] There we write down the relevant terms 
for the volume corrections to the magnetic radius. 
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which produces the total charge. In infinite volume, we can boost this result to an inertial 
frame where the nucleon moves with momentum p = Mn v. By Galilean invariance, we 
expect the current to be j = Qn^v. Using Eq. ((Tj) to calculate the current in this frame, 
however, yields 

(N(p)\J\N(p)) = [Q N - Q N (L)}ev, (10) 

where Qn(L) = LCi(L)/3. The electromagnetic current has been additively renormalized 
by the operator with coefficient C±(L). The theory described by Eq. ([7]) is not Galilean 
invariant. This is the non-relativistic analo g of Lorentz symmetry violation on a torus that 
allows current renormalization described in [28| . The sign in Eq. (J7J) [and consequently that 
in Eq. ({TO]) ] anticipates that the infinite volume current is screened in finite volume. 

We can further apply the effective Lagrangian in Eq. ([7]) to polarized matrix elements. 
For a momentum transfer of q = ^-n, we find 

(N(q) I \J k \N(0) I) - (N(q) | \J k \N(0) |) = \ji N {L) + <W /%(£)] > ( U ) 

M N 

which is sensitive to the finite volume magnetic moment and the additional zero-mode cou- 
pling ~p N (L), which is given by ~p N (L) = eL C2(L)/3. To determine the volume dependent 
couplings, we must match finite volume chiral perturbation theory calculations onto the 
single nucleon effective action in Eq. ([TJ. We now turn to these calculations. 

III. HEAVY BARYON CALCULATION 
A. Chiral Perturbation Theory 

To determine volume corrections to photon-nucleon couplings, we utilize heavy baryon 
chiral perturbation theory. The virtual pion cloud of the nucleon deforms in finite volume 
because pions are the lightest hadrons, and easily propagate to the boundary. Their small 
masses arise because pions are the pseudo-Goldstone bosons arising from spontaneous chiral 
symmetry breaking. For the case of two flavors, the symmetry breaking pattern is from 
SU (2) l <8> SU (2) r — > SU(2)v, and the Goldstone manifold is parametrized by the coset field 
E = £ 2 = exp(2z0//), where the SU(2) matrix <ft contains the pions 



(12) 



At leading order, 0(e 2 ), in a small momentum, p, and pion mass, m n expansion, where 
e r*j p ~ m^, the theory of pions is described by the Lagrangian 

C = ^ [tr(D M E^) + ml tr(S + £t)] . (13) 

o 

Here we work in the isospin limit, and at tree level is the physical pion mass. The 
dimensionful parameter / is the pion decay constant, / = 132 MeV. Electromagnetism has 
been gauged into Eq. (fT3|) via the covariant derivative = + ieA^lQ, ], where Q is the 
quark electric charge matrix 

'2 r 



Q = diag . (14) 
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To include baryons we use the isospinor N from above, 



"=Q. (15) 

as well as the heavy baryon formulation which allows us a consistent expansion in powers of 
residual momentum, ~ e, where we decompose an arbitrary nucleon momentum P into 
P M = M N v^ + k^. To leading order, 0(e), the heavy baryon chiral Lagrangian is 

C = N(iv-D + 2g A S -A)N, (16) 

where A^ and are axial-vector and vector pion fields defined by 

(17) 
(18) 

gauge covariant derivative D^, 

D,N = (d, + V, + ieQ N A„) N. (19) 

The nucleon charge matrix Qn is given by Qn = Q + tr Q. 

The nearest baryon resonances in the spectrum of QCD are the deltas. Phenomenolog- 
ically we know the nucleon-delta axial couplings are not small. Furthermore the nucleon- 
delta mass splitting A ~ 300 MeV is not much greater than the pion mass. In fact on current 
lattices, the lightest lattice pion masses are the same size as A. Hence for a consistent 
power counting, we must treat A ~ e, and retain deltas in the theory as explicit degrees of 
freedom. The deltas are described by the flavor tensor field T M , which is heavy baryon Rarita- 
Schwinger spinor. The leading order, 0(e), Lagrangian for the deltas and their interactions 
with nucleons is given by 

C = (iv ■ D — A) T M + g AN (TM M iV + AU%) + 2g AA T»S ■ AT^. (20) 



the latter appears in the chirally and electromagnetically 
which acts on the nucleon isospinor as 



B. Magnetic Form Factor 

Having spelled out the interactions of nucleons with pions and deltas, we detail the 
calculation of the magnetic form factor in infinite volume. Consider the polarized nucleon 
matrix elements that define the Pauli form factor 

(N(q) i \J k \N(0) |) - (N(q) | \J k \N(0) |) = ^F 2 (g 2 ). (21) 

M N 

The magnetic moment is given by the form factor at zero momentum transfer 

Hn = F 2 (0). (22) 

In the heavy baryon theory, contributions to F 2 (q 2 ) arise from local photon- nucleon interac- 
tions, as well as from virtual pion loop contributions generated by leading-order Lagrangian. 
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FIG. 1: Loop diagrams contributing to the nucleon magnetic moments at O(e). The photon is 
depicted as a wiggly line, mesons are denoted by a dashed line, and a thin solid line denotes a 
nucleon, while the double line denotes a delta. 

The leading contribution to magnetic moments come from the local interactions at O(e ). 
These are contained in the current operator 

J" = {voNiS", S»}N + fnN[S^ S»]t 3 N) , (23) 

M N 

where r 3 is an isospin matrix. The two terms above contribute to the isoscalar and isovector 
moments, respectively. For simplicity we shall address the isovector magnetic moment and 
isovector Pauli form factor. These are defined simply as proton minus neutron differences 
in these respective quantities. 

Higher order corrections to the Pauli form factor come from loop diagrams generated 
from the leading order Lagrangian in Eqs. (|T6|) . and (|20|) . Specifically the one loop graphs 
for the Pauli form factor are depicted in Figure [TJ These graphs make 0(e) contributions to 
magnetic moments. Evaluation of these diagrams, combined with the tree-level result yields 
the isovector Pauli form factor F%\q 2 ) j36|, 37, 38], 



ifV) = 2l*i ~ 4^7? fdxm^P^q 2 ) - 9 -f^f C dx F[m 7T P 7T (x, q 2 ), A], (24) 

Jo 9tt / Jo 



where 



P n (x,q 2 ) = Jl-x(l-x)^, (25) 



encodes the momentum transfer dependence, and the non-analytic function F(m, 5) is given 
by 

FK ^_Mg +V " log ( ^S±| ). (26) 

In order to arrive at this answer we have renormalized the tree-level coupling fij, so that at 
one-loop order the fij in Eq. fl24l) is the chiral limit value. 



C. Finite Volume Form Factor 

We now consider the finite volume modifications to the isovector Pauli form factor. We 
delay the determination of volume corrections to the magnetic moment until Sect. IIVI To 
calculate volume corrections, we use the same setup as described in Sect. [TTl namely we take 
a finite space of volume L 3 with an infinite time extent. The pions, nucleons, and deltas, 
moreover, have periodic boundary conditions which stem from the periodicity of quarks in 
the lattice action, and the fact that hadron fields are point-like objects in the effective theory. 
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To calculate observables in the finite volume theory we use the Lagrangian in Eqs. (fT3l . 
(fIE|) . and (120]) . The action, however, is given by the volume intergral of the Lagrangian 
with periodic boundary conditions on all fields enforced. Thus the same Feynman diagrams 
are generated in the finite volume theory as in infinite volume. The difference is the spatial 
momentum quantization of internal and external lines. We must further choose m^L 3> 1 
so that the zero mode of the pion field does not become strongly coupled 0, (3|. With this 
assumption on the lattice size L, the e power counting we have employed in infinite volume 
carries over to finite volume 

The volume correction to the nucleon current matrix element, 5(N\J fl \N)L, can be de- 
termined by a trivial matching condition, namely 

(N\J»\N) L = (JV|J"|iV)oc + 5(N\J»\N) L , (27) 

where (N\ J^liV)^ is the finite volume current matrix element, and (N\ J^lN)^ is the infinite 
volume result, e.g. the isovector magnetic contribution shown in Eq. (12~21) . We have temporar- 
ily suppressed the momentum transfer argument and denote L-dependence with subscripts 
rather than parenthetically. Because the finite and infinite volume theories share exactly 
the same ultraviolet divergences, the matching condition above ensures that S(N\ J^\N)l is 
indeed the infrared effect. 

Carrying out the matching in Eq. (l2Tj) on the spatial isovector current matrix element, 
we find 

u[S J ,S-q]u, (28) 
where we have defined the function 

f ij (q 2 ,q,5) = I dx{2/g 2 [a;q,m^(x,g 2 ),5]+g^ /2 [xg,m^(a;,g 2 ),5]}, (29) 

J 

in terms of lp(0, m, 5), and l*g(0, m, 5) which are finite difference functions appearing in the 
Appendix. The u, and u appearing in Eq. (1251) are Pauli spinors, and P n (x,q 2 ) is given in 
Eq. ([25]). 

From Eq. (125]) . we cannot identify finite volume corrections to the Pauli form factor. 
Indeed the form factor decomposition of nucleon current matrix elements is tied to Lorentz 
and gauge invariance. The former does not hold on a torus, while the latter has a completely 
different form. The integrand of Eq. (12"§1) is clearly an even function of each component of q. 
More generally it is a cubic invariant function of the spatial momentum, and so depends on 
q 2 , as well as QiQiQiQi: an d a tower of other cubic invariant combinations. As there is no 
simple way to denote this, we have chosen simply to write / lJ = f^(q 2 , q, 5). This function 
should be considered as a finite volume generalization of a form factor. 

With Eq. (125)) . we have determined volume corrections to nucleon current matrix el- 
ements. This result can be directly utilized to remove volume dependence from lattice 
calculations of the nucleon current provided one is in the range of applicability of heavy 
baryon chiral perturbation theory. The expression can be simplified, moreover, depending 
on the actual lattice kinematics used in measuring the three-point function on the lattice. 
For example, to be sensitive to the magnetic coupling in Eq. pip , the momentum transfer 
q must have at least one component transverse to the plane of spin polarization. We shall 
thus assume that q has no component parallel to the spin polarization. Additionally from 
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S(N{q)\JUN(0))i 



P 



Eq. (TTTj) . we see that the spatial component of the current measured must be transverse to 
the spin polarization. We can thus take q and J to be mutually orthogonal in the plane 
transverse to the polarization. This choice is sensitive to the Pauli form factor and simplifies 
the volume dependence. Because q is orthogonal to J, the second term in the integrand of 
f l i(q 2 , q, 5) is zero since i labels the index of the current in Eq. 0281) . 

For the sake of concreteness take q = (q, 0, 0) and consider the y-component of the 
current. Then we have the volume effect that we loosely denote by dF^iq, L) for lack of a 
better symbol, 

5(N(q) | \y ■ J lso \N(0) |) L - S(N(q) | \y ■ J lso \N(0) T)z 

dx j^A^N*, m n P 7I (x, q 2 ), 0] + ^g\ N ll%[xqx, m n P 7I (x, q 2 ), A] j . 

(30) 

If other kinematics are chosen, one must return to Eq. ( 1281) and evaluate accordingly. 



SF^(q,L) = 




IV. MATCHING AND DISCUSSION OF RESULTS 

Having deduced the finite volume modification to nucleon current matrix element, we 
can now analyze the volume dependence and make contact with the single particle action 
in Eq. ([7]). Connection to the single particle action with spin-dependent higher derivative 
terms is carried out in Appendix [Bl 



A. Matching 

To perform the matching, we return to the volume corrections under the general kine- 
matics, Eq. ( 1281) . and series expand in the momentum transfer. To find the leading term, 
we need to evaluate the finite volume form factor f l i(q 2 , q, 5) at zero momentum, 

r(0,0,5) = ^J 5/2 (m.,5), (31) 

where J/3(m,5) is given in Appendix El After some algebra, the leading term in the series 
expansion of the matrix element can be written as 

5(N(q)\JUN(0)) L = ^-u[S\S-q]u6F^(0,L)+O(q 3 ), (32) 

M N 

where F^iO, L) is identical to the function in Eq. (130]) evaluated at q = 0. Comparing with 
Eq. ([7]), we thus find 

^\L) = F 2 (/) (0) + 5F 2 (/) (0,L) (33) 
C 2 (L) = 0. (34) 

The term SF^iO, L) is identified as the finite volume correction to the isovector magnetic 
moment, 5//^(L), and (after algebra) is identical to that determined in j^]. Furthermore 
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the additional magnetic-like zero mode interaction in Eq. (|7j) vanishes to this order. This is 
true of both the isovector and isoscalar combinations. The remaining coefficient in Eq. (JTj), 
C\(L), must be determined by computing the charge form factor in finite volume. The first 
non- vanishing contributions are from recoil order (~ l/M/y) terms in they heavy baryon 
theory. 

The finite volume correction to the current matrix element in Eq. (I28p is determined at 
some general value q = 2nn/L. We must investigate the conditions under which the series 
expansion in quantized q leading to Eq. fl32|) is justified. For simplicity of notation, consider 
the momentum to lie entirely in one direction. A derivative expansion of the matrix element 
is justified as follows. For the n-th mode, q n = 2im/L, and for large enough n, the relative 
difference Sq, 

x qn+i ~ q n 1 / _v. 

Sq = = -, (35) 

q n n 

approaches zero. To rigorously series expand, we must take differences of the nucleon current 
matrix element between adjacent modes for large mode number. This is the procedure 
implicitly needed to arrive at Eq. fl32|) . Such a procedure is lacking when one is restricted 



to just the lowest available momentum mode, q = 2ir/L, see 29|. We shall investigate this 
below numerically. 

We now estimate when Eq. (131Z1) gives a reasonable approximation to the finite volume 
effect. For the effective theory to give a reasonable description of the magnetic moment as 
opposed to the form factor, we must have the higher order terms, 0(q 3 ), under control. For 
~ 20% error in keeping the linear term in Eq. (I32l) . we require 

q2 < I (36) 



Am 2 m 2 L 2 ~ 5 

To have a reasonable approximation to the derivative, say ~ 20%, we must further require 

H=\<\ (37) 

Combining these two restrictions, we find 

m^L > VSirn > 5V5n » 35. (38) 

Strictly speaking, this restriction must be met to determine the volume effect to ~ 20% 
using Eq. (13"2"1) (or using the equivalent expression in j^|) for a lattice determination of the 
current matrix element using the n-th mode. The volume must be large so that there is 
a separation of scales: the mode number must be small enough for the effective theory to 
be applicable, and the mode number must be large enough to allow an approximation to a 
continuous derivative. For volumes this large, however, the finite volume effect can be safely 
neglected altogether. 

The actual lattice determination of the nucleon magnetic moment complicates straight- 
forward application of the finite volume result, 8^{L), derived above. Current lattices 
are obviously too small to allow for the procedure outlined above to extract the moment. 
Instead one is left with quantized modes that are not close enough together so that mode 
differences approximate derivatives. Furthermore the values of the lowest non- vanishing mo- 
mentum transfers are not in the range of applicability of the effective theory, which requires 
q 2 /m% ~ 1. Typically some modeling of the g 2 -dependence of the matrix element is done 
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FIG. 2: Plot of the ratio R(L) in Eq. (j39[) as a function of m n L. 



to extract the moment. It is unclear quantitatively what the volume correction to such 
model-dependent extrapolations should be. In order to shed some light on the subject, we 
analyze the finite volume correction to the current matrix elements, and give a qualitative 
picture for why dF^pfi, L) does not provide an estimate of volume effects relevant to typical 
lattice extractions of the magnetic moment. 



B. Discussion 



Let us assume that the lattice practitioner has determined the isovector current utilizing 
the smallest non- vanishing momentum mode, q = 2ir/L. Consider the magnetic moment 
correction as a correction to the finite volume current rather than just narrowly to the 
magnetic moment. 4 For large enough box size, the smallest momentum mode approaches 
zero, and so we expect SF^iO^L) to give a good description of the volume dependence of 
the finite volume form factor evaluated at the lowest mode. We can address this numerically 
by comparing SF^fi, L) to the momentum transfer dependent volume effect. For ease let 
us assume the lattice kinematics are such that we arrive at Eq. fl30|) in finite volume. We 
plot the ratio R(L) defined by 

R(L) - Wm' {) 

as a function of the box size m^L. This is done in Figure [2j Here q is fixed at the first 
non-vanishing momentum, and we use the values g& = 1.25, |gA7v| = 1-5, M N = 0.94 GeV, 
and m^/A = 0.48. From the figure, we see that <5F 2 (/) (0,L) gives a rough estimate of 
the volume effect for the matrix element under the simplifying kinematical choices. The 
agreement improves as the box size L increases. The ratio R(L) tends to one asymptotically 



4 We shall often refer to current matrix elements and form factors interchangeably. In the present case, 
the F2(q 2 ) form factor is just the current matrix element scaled by q — 2n/L. The same is true of the 
volume effects due to the simplifying kinematics. 
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but only very slowly. Numerically one can see this asymptotic approach set in at m^L ~ 90. 
As explained in 23], the current matrix element depends on two combinations involving 
the momentum transfer: q 2 /m 2 , and q 2 L 2 . For the lowest non- vanishing mode, the latter 
combination is constant and a series expansion is poorly convergent. 

An expansion in q 2 /m 2 for large L is well behaved, and we can obtain a different ap- 
proximation to SF^ (q,L) in the asymptotic regime. It is best to illustrate the issue by 
considering a similar but simpler example. Notice that for large m^L, 

P^q 2 ) L -^l. (40) 
Thus a typical finite volume contribution of the form 

J dxIpi^^x.m^P^x.q 2 ),/^ — > J dx Ip i^-^-x, m^, , (41) 

where we have specified the n-th mode in order to track the problematic term. In this limit, 
the x-integral can now be performed because the only x-dependence enters in one of the 
elliptic-theta functions $3, namely as 

-1 

dx$3(iinx, e~ T ) = 1. (42) 





Thus to take the large L limit of Ip, we merely drop the elliptic-theta function in the direction 
corresponding to the momentum transfer. The resultant volume correction is effectively two 
dimensional. This procedure is only valid when the momentum transfer is aligned with one 
of the spatial axes. The generally oriented case is considerably more complicated but the 
x-integration can still be performed. Expanding in q 2 /m 2 but not in q 2 L 2 , we arrive at 

6F 2 {I) (0, L) = ~f j dx ^g 2 A I 2 %[xqx, 0] + ^g 2 AN l!%[xqx, m n , A] J , (43) 

where we have explained above how the x-integrals are to be performed analytically. From 
this correction, we can form the ratio R(L) = 5F^ ; I \o,L)/5F2 f\q,L). This ratio has the 
same behavior as R(L) shown in Figure [2J The only difference is that R(L) is ~ 10% 
larger. The careful expansion leading to 5F^ ; (0, L) thus actually does worse to describe 
the finite volume matrix element than the naive expansion, 5F 2 (0, L). The asymptotic 
approach to unity also sets in for R(L) at m n L w 90. Of course at these volumes any effect 
is negligible. To compare these approximations to the finite volume matrix element, we 
plot each individually as a function of m n L in Figure [31 For all practical purposes, beyond 
m^L m 5 both approximations give reasonable values for the finite volume correction to the 
current matrix element. 

If the simplifying choice of lattice kinematics is not made, the correction dF^fi, L) 
ceases to give a good estimate for the volume corrections to the current. We can investigate 
this by choosing the momentum transfer to be q' = ^-(1, 1, 0) and returning to Eq. ( f28i) to 

determine the volume effect. This effect we denote by 8F2 (q 1 , L) which has the form 

6F 2 {I \q', L) = ^ \s(N(q') | \y ■ J lso \N(0) [) L - S(N(q') T \y ■ Ji so \N(0) T)l 
Ztti L 

3MN Kl [f 2 W\ 4, 0) - / 21 (g' 2 , q\ 0)] + 2 -gl N [f\q'\ q', A) - f 2 \q' 2 , q> 
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FIG. 3: Plot of the finite volume modifications at q = ^-(1, 0, 0). The finite volume matrix element 
5F 2 I \q, L), the naive guess SF^fi, L), and the approximation SF^fi, L) are plotted as a function 
of m^L. 

Notice we have scaled the matrix elements by a factor of 2it/L, not by \q'\. This is because 
the overall momentum factor in the matrix element is 2ir/L by virtue of Eqs. (|21|) and (|28l) . 
Expanding in large m w L but fixed q' 2 L 2 , we can derive the asymptotics of 5F 2 (?', L), which 
we denote by the function 5F 2 (ft, L) and is given by 

(45) 

The resulting ^-integrals in SF^ffl, L) involve products of two x-dependent Jacobi elliptic- 
theta functions. These can be evaluated analytically as shown in Appendix [A] The notation 
0' reflects that the asymptotic form still depends on the mode numbers (in this case n x = 
n y = 1). We can now compare finite volume effects in the current matrix element for the 
momentum q' = ^(1,1,0). The exact correction 8F 2 (q', L), naive expansion 5F% [ (0, L), 
and asymptotic correction 5F 2 (0',L) are each plotted as a function of m^L in Figure HI 
From the figure we see that naive guess 8F 2 (0, L) ceases to give a good description of the 
volume effect. The asymptotic formula &F 2 (0', L) approaches 5F 2 (q',L) for reasonably 
large values of m n L as it should. There is a crossover for m n L rs 3 which is why the curves 
start out in near agreement on the plot. 

Now we remind the reader that 5F 2 (0, L) was derived as the finite volume modification 
to the magnetic moment, which we denote 8fJjp{L). Above we have been considering this 
term as a correction to the finite volume form factors. The question remains: if we use lattice 
data for the current matrix element, to what extent is the finite volume correction to the 
magnetic moment described by 8^{L)1 Based on our discussion leading up to Eq. (|38j) . 
we expect a poor description of the volume effect with current lattice sizes. We can only 
describe the behavior qualitatively because it depends on precisely how the magnetic moment 
is extracted from the lattice data. 



5F 2 [I) (0',L) = 5F 2 [I \q',L) 
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FIG. 4: Plot of the finite volume modifications at the momentum transfer q' = ^-(1, 1, 0). The finite 
volume form factor SF^p (q', L), the naive expansion 8F2 (0,1/), and the asymptotic approximation 
SFP(0',L), are plotted as a function of m n L. 

Generally speaking what happens in any procedure (including modeling the momentum 
transfer dependence) is some type of weighting over the values of the current matrix element 
for the smallest momentum modes. The simplest way to do this is to take the slope of the 
current matrix elements for the two smallest non- vanishing momenta. By Eq. (I21~|) . the slope 
should be the magnetic moment. Now consider the volume corrections to the current matrix 
element at each of the sampled momenta. For nearby modes, the volume effect has a good 
degree of cancellation and is hence reduced. Compare this to 8F^ T (0, L) which provides an 
otherwise decent estimate of either of the individual values of the current matrix element, and 
consequently a poor estimate of their difference. We expect this argument to qualitatively 
generalize when one considers an arbitrary weighting of contributions from small momentum 
modes. This constitutes a rough description of how one models the momentum transfer 
dependence of form factors. 

The safest route in dealing with volume corrections is to deduce them for the momentum 
modes at which the current has been calculated. Use the effective theory, if possible, to 
isolate the infinite volume values of the form factor at these data points, then perform a 
momentum extrapolation to extract the infinite volume magnetic moment. Carrying out 
the procedure in the opposite order most likely leads to specious volume dependence. 

C. Form Factor Puzzle 

Finally we explore a surprising consequence of our finite volume expressions for current 
matrix elements. We show that values of current matrix elements calculated at two differ- 
ent momenta sharing the same value of momentum transfer squared generally differ. This 
puzzling feature is a manifestation of the lack of rotational invariance in lattice data. 

We consider volume corrections to the isovector Pauli form factor and thus return to the 
general expression in Eq. (1281) . Above we have determined the effect of the finite volume at 
q' = ^-(1, 1, 0). Now we additionally determine the effect at the momentum q" = ^-(1, 0, 1). 
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FIG. 5: Plot of the form factor difference AK 



GO 



, _ F 2 (/) (q' 2 ) -F^(q" 2 ), with q' 2 = q" 2 = 8vr 2 /L 2 , 
as a function of m^L. Rotational invariance at infinite volume requires this difference to vanish. 



-F^\q" 2 ), with q 12 



Of course q = q , and in infinite volume the form factor must be the same by rotational 
invariance. For both momenta, we utilize the ^-component of the current. 5 To compare the 
value of the form factor at these differing momenta, we form the difference AF^ given by 



AF 



CO 



LM 



N 



2m 
LM N 



2m 



(N(q') I \y ■ J lso \N(0) j) L - (N(q') | \y • Ji so \N(0) T)l 
(N(q") I \y ■ Ji so \N{0) l) L - (N(q") | \y ■ J Iso \N(0) T)l 



(46) 



Accordingly the infinite volume pieces of the current matrix elements cancel out of AF 2 



Thus AF 2 (/) gives a measure of how much calculations of the form factor at q = ^jt & n d 
differ. Using the expression for the volume effect, Eq. (1281) . as well as expressions 



,//2 _ 8tt^ 



from Appendix |AJ we determine AF^ to be 



AF 



6M, 



N 



P 



&W) + gSLv/(A) 



(47) 



where the function /(A) is given by 



x [^(ttx, e- r )^ 3 (vrx, e^^O, e~ T ) - ^(ttx, e^f^C), e" T ) - tf' 3 '(0, e^s^x, e^) 2 ] 

(48 



5 A third choice of momentum, q 1 " — ^(0, 1, 1), additionally has the same momentum transfer squared. 
Using the y-component of the current in Eq. (|2ip , this momentum yields zero for spin-polarized matrix 
elements in infinite volume. Any signal from q'" on the lattice would be purely a volume effect. 
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We plot AF 2 in Figure O and indeed find that the finite volume form factor does not 
respect rotational invariance. From the Figure, moreover, we see that for asymptotically 
large volumes rotational invariance is restored. This can easily be demonstrated analytically 
from Eq. (jUJ). This is a general feature of finite volume form factors. One can easily verify 
the lack of rotational invariance for the isoscalar combination as well as charge form factor, 
for example. 

The fact that different momenta (which share the same magnitude) yield differing current 
matrix elements is already evident from the general structure of the finite volume effective 
action, Eq. (J7J) [as well as Eq. (IBlj) in Appendix [B]. The new couplings allowed on a torus 
break rotational invariance. While these couplings vanish at one-loop order in heavy baryon 
%PT, there are indeed non- vanishing couplings at higher orders in the derivative expansion. 
These terms are not suppressed, however, because all terms in the derivative expansion in 
powers of q 2 L 2 are order one. 

In confronting actual lattice data for current matrix elements, we cannot attribute all 
differences at identical momentum transfer squared to volume effects. The discretization 
also breaks rotational invariance. Thus at infinite volume but finite lattice spacing, matrix 
elements are only hypercubic invariant functions of the momentum. Generally differences 
seen in lattice form factor data arise from the violation of rotational invariance from both 
volume and discretization effects. Physically one expects the dominant source of difference 
to be from the volume (discretization) at small (large) momentum transfer. 



V. SUMMARY 

Above we have analyzed the finite volume current matrix elements of the nucleon. Focus- 
ing specifically on the isovector contribution, we have deduced in Eq. fl28l) the finite volume 
correction to lattice calculations of the electromagnetic three-point function. This finite 
volume result does not have an interpretation as a correction to the Pauli form factor. A 
decomposition of matrix elements in terms of Dirac and Pauli form factors makes use of 
Lorentz and gauge invariance. The former is not relevant for a torus, while the latter has 
new features in finite volume due to the special gauge transformation of the zero-mode. 

For the smallest available momentum transfer on the lattice, q = 2ir/L, we have shown 
numerically that one cannot series expand finite volume results in powers of the photon 
momentum. This result applies generally, not just to current matrix elements. Consider 
a simple Feynman diagram with momentum insertion in finite volume. Using a Feynman 
parameter to combine denominators, the typical contribution has the form 



Jo „ L 



An (n + xm) + (m^-L) I 1 — x(l — x. 

mi 



2 \ i -P 



(49) 



where q = 2irm/L is the inserted momentum. In infinite volume, multipole expansions 
are well defined in terms of the expansion parameter A = q 2 /m 2 . In finite volume, this 
expansion remains well behaved, A oc (m^V)~ 2 for large L. There are additional finite volume 
terms, however, that depend on the combination q 2 L 2 = 47r 2 m 2 . Such dependence cannot 
be expanded out even in large volumes. Thus working at only the smallest momentum 
transfer, we cannot deduce volume corrections to multipole moments, cf Figure [2j The 
approximate agreement shown in the figure results from simplifying assumptions about the 
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lattice kinematics used to measure the current. The agreement worsens when such choices 
are not made, see Figure HI 

On the other hand, for large enough momenta the relative spacing between modes be- 
comes small enough to approximate a derivative. In this regime, one can series expand finite 
volume modifications as has been done in the literature. The resulting frequency expan- 
sion is not a multipole expansion, however, but is described by a theory of Wilson lines, 
as in Eq. ([?[). New couplings are allowed and lead to current screening, for example. For 
the isovector magnetic moment of the nucleon, we showed that a new zero-mode coupling 
allowed by symmetry vanishes at next-to-leading order in heavy baryon chiral perturbation 
theory. Consequently we were able to recover previously derived results [§]. In order for 
such results to be valid, however, one requires momenta small enough for the efficacy of a 
low-energy effective theory, and at the same time large enough in order to treat the momen- 
tum as continuous. Satisfying such restrictions is considerably beyond the reach of current 
computing power. In practice, one must retain the momentum transfer dependence of finite 
volume matrix elements to account for the volume effect in lattice correlation functions. 
Finally we showed that such finite volume corrections to form factors can lead to surprising 
effects. Specifically current matrix elements evaluated at identical q 2 , but differing q, need 
not be the same. 
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APPENDIX A: GLOSSARY OF FINITE VOLUME FUNCTIONS 

Above we have determined the finite volume modifications to single nucleon current 
matrix elements. In this Appendix, we give explicit formulae for the finite volume difference 



functions used in the main text. We use similar notation for these functions as [39 . 140 
where further discussion can be found. 

In evaluating a Feynman diagram in finite volume, the loops contain a sum over the 
allowed Fourier modes in a periodic box. The difference of this sum and the infinite volume 
integral is the finite volume effect. As is customary, we treat the length of the time direction 
as infinite. In a heavy fermion formulation, all finite volume differences with momentum 
insertion can be cast in terms of the function 1^ %3 (6, m, 5), defined by 




Ip 1 (0, m, 5) - ^ . .. ^ {q + y + m2 + 25X + X2]l3 

dq q n ■ ■ ■ q l 



(2vr) 3 [(g + 0) 2 + m 2 + 2<5A + A 



21/3 



(Al) 



where n sums over triplets of integers, and the momentum modes satisfy the periodic box 
quantization condition, q = 2irn/L. While general expressions for the exponentially con- 
vergent form of l 1 } lj (6,m,5) exist, we merely cite the required cases for our work. These 
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arc 



1/3(0, m, 5) 



4(0,171,5) 



Ip(0,m,6) 



^T((3) J 



dTT l^ e -(- 2 -^)^/4r Erfc 



5L 



209 



m, 5) - 9^I p (d, m, 5), and 



(A2) 
(A3) 



1 



4(/3-2)(/3- I)d6^d0 
1 



203-1) 



2(^-1) 



+ 9 h 9 t2 I f3 (0,m,5), 
(A4) 



where ^3(2, g) is the Jacobi elliptic-theta function of the third kind, and Erfc(z) is the 
complement of the standard error function. The finite volume physics arising from the mass 
splitting parameter 5 is well explicated in (Io| . 

In expanding matrix elements about zero external momentum, one also has use for the 
functions Jp(m,S), and Kp(m,5) given by 



Jp(m,6) 



1 



& 2 



A((5 - l)((3 - 2) 88 2 

1 d 2 
' 4((3 - 1)((3 - 2) 85 2 



i>_ 2 (0,m,<J) + 1 + 



Jp- 2 (0,m,5) + 1 + 



(5 d 
(3 -IdS 

5 d 
(3-186 



1/3-1(0,771,5) - m 2 1^(0,171, 5), 

(A5) 

7/3-1(0,771,5) - m 2 J(3(0,m,5). 



(A6) 



Lastly of use are Feynman parameter integrals of Jacobi elliptic-theta functions. 

•1 







dxf) 2 


>x, g) = 


1, 


(A7) 


dxf) 2 


>x,g) 2 = 




(A8) 


dx tfg 


[irx,q) 2 = 


= -^(o,g 2 ), 


(A9) 


da #3 


{TTX,q)^ 


;7rx,g)=^(0,g 2 ). 


(A10) 



Here the primes denote derivatives with respect to the first argument. Double primes can 
be traded in for single derivatives with respect to the second argument because the "&z(z, q) 
satisfy a diffusion equation. 



APPENDIX B: MAGNETIC RADIUS 

In this Appendix, we present the extension of the finite volume effective action in Eq. ([7]) 
to higher order in the derivative expansion. This will allow us to consider the finite volume 
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corrections to the magnetic radius. For simplicity, we ignore the charge radius and focus 
just on the spin-dependent part of the single particle effective action. 

Using the cubic symmetry of the torus, it is easy to see that there are no spin- dependent 
terms with two derivatives. Including all spin-dependent terms with three derivatives and 
at most one insertion of VV„- , we arrive at the nucleon effective action 



N 



C 3 {L)a ■ d 2 B - C 4 (L)tr ■ V 2 B - C 5 {L)(T ■ (V x V 2 W 



-C 6 (L) ^(TiViViBi - C 7 {L) J^etffcOiViViVjWj 



(-) 

k 



N. 



(Bl) 



Here the couplings C^(L), C^(L), Cq(L), and C 7 (L) run to zero as the volume goes to 
infinity. This reflects that the corresponding operators are forbidden by infinite volume 
gauge invariance, and Lorentz invariance. The remaining coupling C 3 (L), satisfies 



C 3 (L) = C 3 + SC 3 (L), 
where C 3 is the infinite volume coupling given by 



C 3 = 6 < r 2 M > 



with < r 2 M 



(B2) 



(B3) 
The 



> as the mean-square magnetic radius which we define by —^-^^(q 2 
finite volume effect 6SC 3 (L) is thus the correction to the magnetic radius, however, there 
are four additional operators that also make contributions to spin polarized current matrix 
elements in finite volume. 

To determine the couplings in Eq. ( 1B1I) . we turn to heavy baryon chiral perturabtion 
theory. Specifically we need the result derived above for spin polarized current matrix 
elements, Eq. (1281) . We expand this result to third order in the photon frequency and match 
onto Eq. (IBlj) . This will determine the isovector part of the coupling constants. We find 



5C 3 (L) 

C 4 (L) 
C 5 (L) 



6/ 2 
1 



g\J 7/2 {Tn^, 0) + -J 7/2 (m n , A) 



9a (7^9/2^, 0) - 5J 7/2 (m^, 0)) + -g 2 AN (7K 9/2 (m n , A) - ^Jj^m^, A)) 



3/ 2 

C 6 (L) = C 7 (L) 



0. 



(B4) 



The finite volume functions Jp(m,5) and Kj3(m,5) are defined in Appendix El The fact 
that C 3 (L) and C±(L) are non- vanishing to this order means that it is not possible to define 
the magnetic radius at finite volume. The remaining coefficients require further terms in the 
chiral expansion to be non-vanishing. The non-relativistic nature of the nucleon in heavy 
baryon chiral perturbation theory is likely the culprit for the vanishing of these terms at 
next-to-leading order. In the meson sector, by contrast, derivative couplings lead to a myriad 
of new finite volume terms, see 12911. 
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